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Abstract 

One of the less well understood ambiguities of quantization is emphasized to re¬ 
sult from the presence of higher-order time derivatives in the Lagrangians resulting in 
multiple-valued Hamiltonians. We explore certain classes of branched Hamiltonians in 
the context of nonlinear autonomous differential equation of Lienard type. Two eligible 
elementary nonlinear models that emerge are shown to admit a feasible quantization 
along these lines. 
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1 Introduction 


The Bohr’s formulation of his model of atom [T] was one of the most important steps towards 
the subsequent birth of full-fledged quantum theory from classical mechanics. His choice of 
benchmark systems was insightful (quantum hydrogen with Coulomb potential V(x) l/|f| 

in non-relativistic Schrodinger equation is one of not too may exactly solvable real-world 
systems) and lucky (indeed, the subsequent inclusion of subtleties led just to some not too 
essential corrections and modihcations of the picture). 

Today, more than 100 years later, the Bohr’s emphasis on connectedness between 
classical and quantum worlds is still alive and inspiring. Although we are currently paying 
attention to multiple much more sophisticated quantum systems, the principle of corre¬ 
spondence and the traditional description of dynamics via certain real and local potentials 
V (x) still guides the majority of phenomenological studies in atomic, molecular and nuclear 
physics. 

In a broader methodical framework, say, of condensed matter physics, the traditional 
approach does not work that well. Incessantly, its applicability is being challenged by the¬ 
oretical innovations. Just to name a few, let us mention the recent growth of popularity 
of the detailed study of certain non-real, complex potentials with real spectra [2] or, most 
recently, of their friendly and tractable non-local generalizations [3]. An emerging interest 
also involves various non-linear forms of the interactions which may be mediated, say, by an 
energy-dependence of the dynamical parameters |1]. Last but not least, let us mention the 
fruitful transfer of the dynamics-determining role from elementary potentials V{x) to the 
other factors like, e.g., coordinate-dependent masses m = m{x) [SIEIE] etc. 

In this context a new class of innovations of the description and simulation of quantum 
dynamics recently emerged in connection with the possible specihc role of the models of 
physical systems in which the Lagrangians possess the time derivatives in excess of quadratic 
powers iiiniiiniiiiiinKia]. The use of these models leads, on both the classical and quantum 
level, to the necessity of a re-evaluation of the dynamical interpretation of the momentum p 
which, in principle, becomes a multiple function of velocity v. In what follows we intend to 
present a few results in this direction. 


2 Branched Hamiltonians 

A typical classical model of the above-mentioned non-quadratic type may be sampled by 
Lagrangian 

L = (u-1)^-H(x). (1) 

Function V (x) represents a convenient local interaction potential while the traditional kinetic- 
energy term is tentatively replaced by a fairly unusual function of “velocity” v. 

This model was recently analyzed in Ref. |T3] where the dehnition of the fractional 
powers of difference v — 1 was adapted to the needs of possible phenomenology. In detail, the 
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{2k + 1)—st root was required real and positive or negative for n > 1 or n < 1, respectively. 
Correspondingly, the quantity v turned out to be a double-valued function of p. 

2.1 The problem of quantization 

In a search for a quantum version or analogue of model CO. both of its latter features 
appear truly challenging. Firstly, the dehnition using multivalued functions (referring to the 
literature let us speak about “branched Hamiltonians”) seems to require a deeper insight in 
its physical nature and connections. 

In this note we decided to consider the subject of the branched Hamiltonians in the 
context of a class of nonlinear autonomous differential equations of Lienard type [I5l [16] since 
these systems are of potential importance not only in optics ini but also in the shallow- 
water-wave studies [18] and in non-Hermitian quantum mechanics [19]. Due to nonlinearity 
one of the main obstacles lies in the absence of superposition principle. The recent progress 
in the direction of seeking analytical and numerical solutions as well as of tracking down new 
mathematical properties has been dramatic, nevertheless. In particular, an inspiring feature 
of non-linear oscillations was found in the dependence of amplitudes on the frequency, etc. 
A useful role in this analysis was played, e.g., by the study of the class of Lienard equations 

X + r{^x)x + q{^x) = Q ( 2 ) 

where overdot indicates a derivative with respect to the time variable and where r(a;) and 
q{x) are two continuously differentiable functions of the spatial coordinate x. 

After a restriction of the model to the specihc form of r(a;) = kx and q{x) = Aa:-1- ^x'^ 
(which are odd functions of x) one obtains 

k‘^ 

X + kxx H- x^ -1- Ax = 0, A > 0 . (3) 
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This equation represents a cubic oscillator subject to a damped nonlinear force as indicated 
by the product kxx. Naturally, the presence of the damping is a challenge whenever one 
tries to contemplate a quantization of the model. 

2.2 Hamiltonians 

The classically damped systems leading to non-Hermitian quantum Hamiltonians need not 
necessarily imply a quantum damping [ 20 ] • Moreover, a backward classical-optics reinter¬ 
pretation of quantum Hamiltonians of this type may exist in photonics playing the role of so 
called lossy amplihers [21]. Last but not least, the similar classical-physics contexts (related, 
e.g., to the study of Bose-Einstein condensates [221123]) re-open the phenomenological appeal 
of classical nonlinearities as sampled by Eq. ([1]). For all of these reasons it is worth checking 
the conversion of the Lagrangians of Eq. ([3]) into Hamiltonians. 

A brief account of some background and references is in order. First of all, let us 
recollect [ 21 ] that the toy-model Lagrangians of the form 
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( 5 ) 


2A:2 ' 3 ' 2k 2k^ 

admit the evaluation of the canonically conjugate momentum, 

27A^ . 2 3A 

Thus, the Hamiltonian can be cast in compact form 


H 


_ ^ 


2-2 



k'^x^ 



2kp 2k^x^p 
~ 27A2 


( 6 ) 


Such a Hamiltonian is of non-standard type. The co-ordinates and potentials are mixed 
so that the expression cannot be split into individual kinetic and potential energy terms. 
However, we can write 


+ (7) 

where 

/(p)=.^(i-g). ,8) 

with uj = \/A. The roles of coordinate and momentum have been transposed implying a 
momentum-dependent system at play. The first term in (0) represents a mixed function 
of both position and momentum variables while the second term is a function of momen¬ 
tum alone. The classical Hamiltonian 77(x,p) is invariant under a joint action of coordinate 
reflection and time reversal transformation. 

Exact quantization of the Hamiltonian can be carried out by going over to the mo¬ 
mentum space with x = ih-^. The Hamiltonian turns out to be of momentum-dependent 
mass type. Adopting a von Roos strategy of quantizing the problem by considering a general 
symmetric ordering |25], the underlying Schrodinger equation can be formulated and solved, 
in principle at least. 


3 Two illustrative models of dynamics 

In practice, the price to pay for having the model quantized lies in the clarihcation of the role 
of singularities in the eigenfunction. In this sense the road to quantization is, expectedly, not 
free of difficulties. One can, however, feel encouraged by the work by Chithiika Ruby et al 
[21] who showed that the energy spectrum and normalized solutions could still be obtained 
for a class of Hamiltonians that are nonsymmetric and non-Hermitian. 

Needless to add that the non-Hermitian nature of quantum Hamiltonian may bring 
about a number of unpleasant consequences as, for example, the emergence of the exceptional 

















points [26] or the breakdown of the adiabatic theorem [271 128]. At the same time, the 
acceptance of the anomaly may prove innovative, e.g., by giving a new physical interpretation 
to wave packets [29] or to pseudospectra [HHl 1^ . 

3.1 An amendment of the Lagrangian 

In a way inspired by Eq. ([1]) let us consider a higher-power Lagrangian 

L = C{v + f{x)f^ -6, C= (5)1^, 5>0. (9) 

The main difference from ([T]) lies in our choice of a general function f{x) in place of f{x) = — 1 
in our initial L. The other point is that the inverse exponent with respect to the model of 
Curtright and Zachos na is done for convenience of calculus. Further, we have omitted 
an explicit presence of the potential function assuming that the interaction re-appears in 
a more natural manner via a suitable choice of an auxiliary free parameter 6 and of a 
nontrivial function f{x). As long as our Lagrangian L is of a nonstandard type, we will not 
feel disturbed by the absence of the explicit potential V{x). 

Parameter C is non-negative for 0 < m < | and the canonical momentum is given by 
formula 

dL 2 2 

p =—(5) 1-2™ (n-1-/(x)) (10) 

ov 

which can easily be inverted to yield 

v =-f{x)^. (11) 

The Hamiltonian has the structure 

H±(x,p) = (-p)/(x) - +5. (12) 

2m -|- 1 

For the present purposes it is worthwhile to inquire into the specihc case with m = 0. 


3.2 Type I Hamiltonian at m = 0 

At m = 0 we easily derive the double-valued 

V =v± = -f{x) ± —^ . 

The Hamiltonian branches out to the components 

H±{x,p) = {-p)f{x)T26^/^ + 6. 


(13) 


(14) 
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It is readily noticed that the real or the complex character of H± depends on the sign of the 
momentum p. Once we specify 


„ , A 2 9A" 

/(^) -r 


9A^ 


then under a shift p -^p — 1, H± move over to the corresponding forms 

2 ^ 2(1 


9A^ 

H±{x,p) = ^ 


2kp\ 2 ^ 2kp 2k^x‘^p 


3A 


9A 


3A 


27A2 


(15) 


(16) 


These represent a set of plausible Hamiltonians for the nonlinear system ([3]) possessing 
coordinate-momentum coupling terms. While coincides with the form of hf(a;,p) as fur¬ 
nished in ([6]), in both H± the presence of a linear harmonic oscillator potential is revealed in 
the limit A; —?• 0. The specihc choice of the parameter 6 as given above allows equation ([3]) to 
be converted to a harmonic oscillator form under the nonlocal transformation U = xes 

ra- 

In the classical scenario one needs to restrict p to —oo < p < ^ in order to address 
the physical properties of the system in the real space. However the presence of a branch 
point singularity at p = |^ makes the study of H± quite complicated. At the boundary 
p = 1^, the two Hamiltonians H± coincide. 


3.3 Type II Hamiltonian 

We next focus on a class of models under the influence of the Lagrangian 


L 


1 

s 


^ 2 \ 
-sx -A - u 
6 S J 


-1 


where s is a real parameter. This implies, for the canonical momentum, the result 



1 

s 




(17) 


(18) 


Inversion of (18) leads to 

1 9 3 , 1 

V = -sx^ -h -A ± —= . (19) 

3 s y/sp 

Finally, employing the Legendre transform, the accompanying Hamiltonian to the above 
Lagrangian can be projected, for the two-u branches, as 

H±{x,p) = pu±(p) - L = ^x^p+ -Ap=F 2W- . (20) 

3 s y s 

Eq. (20) describes a set of alternative Hamiltonians for (3). The component iL+ for the case 
s < 0 has earlier been derived in [32] and found to possess an interesting A —?• 0 limit. 
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4 Quantization 

4.1 Schrodinger equation 

Quantization of the branched system ffTBl) can be handled in a typical way as detailed in 
PI, i.e., by adopting a suitable ordering procedure and implementing appropriate boundary 
conditions. In order to carry out quantization for the Type II Hamiltonian we will consider 
here the case of s > 0 and subject it to a perturbative treatment. In this regard we will 
study the spectrum of Hamiltonians entering Schrodinger equation 

H±{x,p) (21) 

where, in units h = 1, 

5/?,n.p)/p=-T:p6yX+2A, (22) 

After the change of variables p = (with special g = 2) and having set pip^'^^p) = r^x(r) 
(with special ■C = 2 + |) one ends up with the following transformed Schrodinger equation 

, e(e + l) , , 36 2 p X 24 , , 12 ^ , 

- X (^^) +- 5 — X{r) + ^r x(r) T^r x{r) = —p'- > x{r) (23) 

where the emergence of the centrifugal term is a byproduct of the change of variables [33] 
and where the quantity resembling the angular momentum is equal to i = 1/2 for our choice 
of exponent g = 2. 

4.2 Perturbation expansions 

Eigenvalue problem fl23]l is well suited for application of Rayleigh-Schrodinger perturbation 
theory [31]. Towards its implementation we abbreviate = gV. Assuming here that 

the value of s is sufficiently large, the auxiliary variable g may be treated as a small parameter. 
Hamiltonians H = as well as the related inhnite sets of wave-function kets |x) = 

and of the bound-state energies E = may be then formally Taylor-expanded, 

CO OO 

lx) = ^ s'” Ixm). E = J29"'E.„. (24) 

m=0 m=0 

Recalling also that the unperturbed Hamiltonian stands for the wave harmonic 

oscillator living on half-line of r G (0, cxo) we arrive at the well known formulae for the 
zero-order perturbation energies Eq = uj (4n -|- 2^ -|- 3) where u = 6s“^. The unperturbed 
eigenfunctions may also be expressed in the compact form 

(^Ixo) = AAr^'''^exp(—a;r^/2) L/"^(a;r) (25) 
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containing a suitable normalization constant M and associated Laguerre polynomials L^n'^ {oj r) 
with, in our present special case, a = I + 1/2 = 1. 

At this moment one would have to return to our hypothesis of smallness of g and 
to prove that the inhnite-series ansatz of eq. fIMl) are all convergent whenever \g\ < R, 
i.e., at some non-vanishing radius of convergence R > 0. Fortunately the proof is an easy 
consequence of the fact that our operator of perturbation V is relatively bounded with respect 
to 

In the final step of the construction of bound states we only have to insert formulae fl2^ 
in Schrodinger equation, deduce the recurrent Rayleigh-Schrodinger perturbation recipe and 
write down the corresponding general formulae that yield, typically, 

= 2n 2 =F ri[n) n = 0,1,... (26) 

where the quantities ?](„) must be evaluated numerically. 

5 Summary 

Although virtually all of the available quantization recipes are ambiguous, one can still 
restrict one’s attention to certain sufficiently elementary classical scenarios and succeed in 
the construction of a mathematically consistent quantum model the practical acceptability 
of which is to be verified via its appropriate experimental realization. 

In such a setting the key motivation of our present study emerged after we imagined 
that the classical double-valued Hamiltonians could offer a unique possibility of moving over 
to the quantum sector thus opening up branched partners in such a setting too. We decided 
to follow this idea and to perform the hrst steps in this direction. 

Our present results should be read as encouraging, confirming the purely technical 
feasibility of the general strategy we described. During its application we revealed that there 
exist several nontrivial examples of the branched Hamiltonians for which the quantization 
appeared not only feasible but also comparatively friendly. 

Beyond this framework, our present approach also seems to open an entirely new 
direction of research in which one might be able to combine the generic numerical techniques 
with the non-numerical insight provided by the mathematics of multivalued functions in 
combination with the analytic tools of perturbation theory. 

In the future research we expect that our present approach could also open a way 
towards the emergence of new qualitative features of branched Hamiltonians. In particular, 
we expect that the new classes of doublets of supersymmetric-partner Hamiltonians could 
emerge here in some exotic manner and in a way which has not yet been followed in the 
literature (cf., e.g., [3^ [36l l38] b 
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